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I. INTRODUCTION 
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The distribution function of local amplitudes, t = \ip(r )\ 2 , 
of single-particle states in disordered conductors is calculated 
exactly on the basis of the supersymmetric cr-model approach 
using a saddle-point solution of its reduced version. Although 
the distribution of relatively small amplitudes can be approx- 
imated by the universal Porter- Thomas formulae known from 
the random matrix theory, the statistics of large t's is strongly 
modified by localization effects. In particular, we find a mul- 
tifractal behavior of eigenstates in 2D conductors which fol- 
lows from the non-integer power-law scaling for the inverse 
participation numbers (IPN) with the size of the system, 
Vt n <x £-(«-!)<!*(«), w here d*(n) =2- f3~ 1 n/(4 : n 2 vD) is a 
function of the index n and disorder. The result is valid for all 
fundamental symmetry classes (unitary, /3 U = 1; orthogonal, 
Po = | , and symplectic, f3 s = 2) . The multifractality is due to 
the existence of pre-localized states which are characterized 
by power-law envelopes of wave functions, |^t(r)| 2 oc r~ 2 ^ , 
[i = fi(t) < 1. The pre-localized states in short quasi-lD wires 
have the power-law tails |V>(a:)| 2 oc x~ 2 , too, although their 
IPN's indicate no fractal behavior. The distribution function 
of the largest-amplitude fluctuations of wave functions in 2D 
and 3D conductors has logarithmically-normal asymptotics. 



Localization of a particle by a random potential has 
been extensively investigated during the past several 
decades It is well known |s|-|7) that, at strong disor- 

der, single-particle wave functions are confined and have 
exponentially decaying tails beyond the scale of the lo- 
calization length L c . At weak disorder, the localization 
length can be very large in ID and 2D conductors, and 
infinite in 3D. A natural question arises: What is the 
behavior of the wave functions at distances smaller than 
the localization length? Despite of its importance, the 
problem of structure of quantum states of weakly disor- 
dered conductors for scales below the length L c has only 
recently started to attract interest p|-|l7|]. 

In particular, one of the issues that has not been ex- 
plored up to now concerns the way the localized states 
develop as a consequence of the increase of disorder in an 
isolated piece of a metal, though a great deal is already 
known about the extended states in it. Some part of 
recent results related to extended (metallic-type) states 
have been obtained by mapping the problem of quan- 
tum mechanics in the classically chaotic systems to the 
Wigncr-Dyson random matrix theory 
the zero-dimensional supermatrix cr-model |7,t 
which are two equivalent ways of describing disordered 
and chaotic systems. 

Both the advantage and disadvantage of such an ap- 
proach comes from the statistical equivalence of eigen- 
states which is usually built into the construction of the 
random matrix substituting the real dynamics. This re- 
veals the set of universalities of the spectra, the level- 
level correlations and the transition matrix elements |l9| ] 
which are similar for a wide variety of objects. For ex- 
ample, the distribution function of local densities of wave 
functions ^(fo)! 2 in a chaotic cavity which one can find 
in such a way is determined only by the fundamental 
symmetry of the system and its volume V ~ L d , but is 
independent of the level of disorder (i.e., of the value of 
a mean free path I) or a physical dimension, d. 

On the other hand, this approach hides individual fea- 
tures of physically different systems and permits to de- 
scribe only metallic-type states which equally test the 
random potential all over the sample. More complex 
states which can distinguish between the ballistic and 
diffusive regimes have to be analyzed beyond the con- 
ventional random matrix theory. Numerical evidence 
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for their existence have been obtained by several groups 
|To| |l2|,0 ■ The goal of the theory to be presented in the 
present paper is to find manifestations of these precursors 
of localization among the wave functions of classically dif- 
fusive conductors {ppl ^> 1, I -C L). That is, we consider 
an isolated piece of a disordered metal with dimensions 
I -C L < L c , assuming that the internal 'conductance' g 
which one would assign to the 'electric circuit' connecting 
the observation point (blown up to the mean free path 
size) with the external surface of the specimen p8| is 
much larger than the conductance quantum, i.e. g 1. 
In the following Sections, we perform a statistical analysis 
of local densities and partly reconstruct spatial structure 
of those rare states which have locally too high ampli- 
tudes (as compared to the average V^ 1 ) to fit into the 
universal random matrix theory description. As can be 
suggested on the basis of the theory below, these states 
are responsible for non-Gaussian tails of distributions of 
fluctuations of local densities of states and conductances 
suggested by Altshuler, Kravtsov and Lerner ||] and are 
generic for the long-living current relaxation discussed in 
Refs. [10. 

Our paper is organized as follows. In Section II, we 
introduce the notion of the eigenstates statistics (II. A), 
discuss the universal distributions of metallic-type states 
(II. B) and, then, sketch the main results of the paper 
focusing our attention at the localization effects (II. C). 
Sections III and IV are devoted to the presentation of our 
theoretical scheme: We derive a reduced supersymmctric 
a- model and show that it has a non-trivial saddle-point. 
The details of the derivation of the saddle-point solutions 
of the reduced cx-model are given separately for each of 
the fundamental symmetry classes (unitary, orthogonal 
and symplectic) in Sections IV.A-C, and the influence of 
fluctuations around the saddle-point is discussed in Sec- 
tion IV. D and the Appendices. The resulting statistics 
of wave functions and the structure of the pre-localized 
states in the conducting regime in ID, 3D and in the most 
interesting case of 2D samples are discussed in Sections 
V, VI and VII, respectively. Section VIII contains a brief 
summary of our results and their discussion. 

II. METALLIC VERSUS PRE-LOCALIZED 
STATES IN THE EIGENSTATES STATISTICS 
(PRELIMINARIES AND RESULTS) 

In this Section, we give mathematical formulation to 
the problem of the eigenvalues statistics in disordered 
systems and consider alternative approaches to its solu- 
tion. That is what is the subsection A about. The next 
part B is devoted to the universal statistics of metallic- 
type of states known in the random matrix theory as 
the Porter-Thomas distribution. The localization effects 
which are beyond the random matrix theory approach are 
discussed in the subsection C, where we give an essence 
of the obtained results. This subsection is written for the 



first reading and can be used as a guide through the rest 
of the text. 



A. Definitions of the eigenstates statistics 

To define the statistics which we shall be studying in 
this paper, we first mention that the properly normal- 
ized eigenstates {ip a } we consider below correspond to a 
quantum particle in a disordered cavity 

{ 5m + u{r) } ^ (r) = £q ^W' ^ r e 5 ) = °' 

where U is a random potential. The local amplitude ip of 
a wave function at some observation point r Q inside the 
sample, i.e., 

t=Wo)\ 2 , (1) 

will be the object of our statistical analysis. In this we 
employ two related quantities: the distribution function 
f(t) of local amplitudes t averaged over disorder, 

f(t) = A{j: a 8(t-\Mro)\ 2 )5(e-e a )) : (2) 

and the set of generalized inverse participation numbers 
(IPN) ||jn| which are the moments of the distribution 
function /, 

t n = A (£ a |V> Q (r )| 2 "<5(6 - e a )) = / °° t n f(t)dt. (3) 

As indicated, () denotes the averaging over random con- 
figurations of a random potential U in the system. In 
Eq. (j|), the sum is over the full set of states {4> a }, V is 
the volume of the system, and A = (vV)~ denotes the 
mean level spacing with v = v(e) the density of states 
per unit volume. Since the distribution function / and 
wave functions {ip a } are normalized, one has the follow- 
ing relations: 

POO 

< = / /(*)* = 1; ti = (iV'al 2 ) = V~ x . 
Jo 

One can also introduce the distribution / s (<7, t) of a local 
spin-density of the wave with a =[,] . It is an important 
quantity for systems with a strong spin-flip scattering. In 
random matrix theories these are known as a symplectic 
ensemble. In such a case, the statistics can be formulated 
in terms of spin-projected eigenstates, e.g. t± — \^i(r )\ 2 
Iplf . The distribution of a total local density (t = t[ 
can be found as the convolution 

f s (t) = f*fa,t-t>)f(i,t>)dt>. (4) 

Historically, the studies of eigenstates in disordered 
conductors started from Wegner's perturbative calcula- 
tions of IPN's Q. Due to the equivalence between the 
descriptions based on the distribution function and the 
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full set of its moments |29|, in most of the later stud- 
ies ^,^,^TJ the eigenstates statistics were reconstructed 
from the set of IPN's. Alternatively, one can start 
from c alculating directly the entire distribution function 
|33], especially regarding the possibility to ap- 
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ply the supersymmetry technique fj. This alternative 
approach has already been used for describing the eigen- 
states statistics over the entire crossover regime from the 
orthogonal to unitary ensembles (low magnetic fields) 
|p6f . More recently, this construction has been advanced 
by developing a reduced cr-model which is applicable to 
closed systems. The reduced cr-model has non-trivial 
saddle-point solutions which enabled us to consider the 
localization effects non-perturbatively |32| . 

The idea to work with the distribution function as a 
whole has also the additional advantage that it makes 
possible to select those rare states which do not fit to 
the universal statistics and study their spatial structure. 
The latter information is implicit in the cross-correlation 
function R(t, r), 

R(t, r) = A (J2 a $ (* - lV^(r )| 2 ) |V^(r + r)\H(e - e a )) . 

(5) 

As it will be clear from the calculations below, rare pre- 
localized states show up as deviations of the function in 
Eq. (0) from universal distributions at the tails where 
t V~ 1 g, so that the combination R(t,r)/ f(t) mimics 
the envelope |V't( r ')| 2 of these states at distances r = |r| > 
I from the top-amplitude (t) position. 



B. Metallic states and universal statistics 

To find the manifestation of the pre-localized states in 
the distribution function /(£), we have, for comparison, 
to give an idea about what would be the form of the 
distribution function if all states were extended. Quali- 
tatively, the extended states test the realization of a ran- 
dom potential equivalently all over the sample and that's 
why their statistics coincides with the Porter-Thomas 
eigenstates statistics renown in the random matrix the- 
ories pTj| , p^ |. The recent studies of the properties of 
the eigenstate of disordered and ballistic chaotic cavi- 
ties (both using the numerical tools (2^j2^] or based on 
the zero-dimensional limit (0D) of the supcrsymmetric 
non-linear cr-model (0,^,^,0) have confirmed such an 
expectation. 

Depending on the fundamental symmetry class, the 
Porter-Thomas distributions can be represented as fol- 
lowing. For the single-particle Hamiltonian describing a 
spin-less particle in the system with a broken (e.g., by 
a magnetic field) time-reversal symmetry (unitary class), 
the distribution function of local amplitudes and the cor- 
responding IPN's have the form 



whereas in the case of a system with the time-reversal 
symmetry (orthogonal ensemble), 



fo(t) 



exp{-^} 



t 



(2«-l)!! 



(J) 



For spin-i particles which undergo a strong spin-orbit in- 
teraction (symplectic ensemble) the Porter- Thomas-type 
of a distribution can be repeated both for the spin- 
projected wave functions and for the total density and 
has the form 



-2Vt 



(8) 



C. Localization effects and eigenstates statistics 
beyond the universal limit 

The universal statistics described by Eqs. (^j^) are 
presented only as a reference point for the subsequent 
analysis. The rare events which cannot be described 
on the basis of the statistical equivalence of eigenstates 
have to be studied using more sophisticated methods, and 
need the extension of the non-linear cr-model beyond the 
OD limit. Details of these calculation are presented in 
Sections IV- VII, whereas in the forthcoming subsection 
we sketch only the basic results. In few words, the uni- 
versal disorder-independent laws work well enough either 
until this disorder is so weak, that the system behaves 
as in the nearly ballistic regime, or at small amplitudes 
t < V ■ s fg. But they are partly broken or, at least, 
modified after the disorder makes the electron motion 
diffusive. In one- and two- dimensional conductors, this 
requires a different statistical treatment of states which 
have too high splashes of a local amplitude, t > V^^/g. 

The method of taking into account all finite (i.e., not 
only small) inhomogeneous variations of the fields used 
in the supersymmetric field theory is based on the exis- 
tence of saddle-point solutions of the non-linear cr-model 
discovered by Muzykantskii and Khmelnitskii jL7|. The 
existence of the saddle-point solution is especially promi- 
nent for a reduced version of cr-model formulated and 
solved in Ref . |3^] . An interesting result of Ref . |32| for 
the quantum diffusion in the dimension d = 2 is the mul- 
tifractality of the states which is in agreement with previ- 
ous numerical simulations (HUHI . Using the same saddle- 
point method as for the unitary (u) symmetry class [p2| , 
we extend the analysis of 2D systems to the other fun- 
damental symmetry classes (orthogonal, o; symplectic, s; 
spin-unitary |34|, su) and arrive, again, at the multifrac- 
tality. The latter is manifested by the following scaling 
of INP's, 



1 



/„(*) = Vexp{-Vt}, t n = n\V~ 



(6) 



Vt n cx L-("-D d *, d*(n) = 2 - f3 = A . (9) 

fts,su 2 

The fractal (or generalized Renyi Q) dimensions, d*(n), 
obey Eq. (||) only for those n's where they are positive 
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and are obtained in the leading order in (2-itvD)^ 1 : where 
D is the classical diffusion coefficient. The sensitivity of 
the derived statistics to boundary conditions, as well as 
the form of the correlation function R(t, r) which we find 
in our calculations, enables us to suggest such a behav- 
ior of 2D multifractal states which is associated with the 
power-law envelopes of their tails, \ip(r + r)| 2 oc (l/r) 2fJ -. 
Being extended from the position of a rare high ampli- 
tude splash |?/;(r )| 2 = t » 1, these tails have exponents 
n(t) individual for each state marked by its own t. 

The behavior of pre-localized states in a quasi-lD wire 
within the localization length scale L c also resembles the 
power-law localization. Independently of t, \ip(x +x)\ 2 oc 
L c x~ 2 . The density of wave functions accumulated by 
these tails is integrable, so that no assertion about frac- 
tality can be made and the inverse participation numbers 
t n tend to take a volume- independent form for n > g. 

The localization effects in 3D conductors are weak, if 
disorder is weak enough to keep the system far away from 
the metal-insulator transition, ppl 3> 1. As a result, the 
statistics of eigenstates in 3D conductor is most similar 
to the universal one: Amplitude t scales with V, and the 
mean free path I appears only as an extra parameter both 
in the distribution function / = f(Vt,ppl), and inverse 
participation numbers t n cx V~ n K(n,pFl)- 

Nevertheless, even then, the statistics of rare events 
shows intriguing deviations from the Porter-Thomas for- 
mulae. Both for 2D and 3D diffusive samples we obtain 
the logarithmically-normal distribution of large local am- 
plitudes of wave functions, 



/(t)«exp{-/3^f ln 2 r}, T 



VtVd 

27T 2 !/D ' 



(10) 



where r/2 = In j and 773 ~ (2£) _1 . Although we study 
an isolated specimen, this result strikingly resembles the 
asymptotics of distributions of local density of states or 
conductance fluctuations found in open systems || . This 
signals about deep physical reasons behind it related, 
probably, to the properties of random walk paths. 



III. EIGENSTATES PROBLEM IN TERMS OF A 
NON-LINEAR SUPERSYMMETRIC a-MODEL 

In the following paragraphs, we formulate the eigen- 
states statistics problem in terms of the non-linear a- 
model. The details of this technique are described in the 
review article [Q, and below we give only a compressed 
extraction from it, keeping similar notations. 

One can try to use the supersymmetry technique as 
soon as a physical quantity of interest is expressed in 
terms of retarded and advanced Green's functions, 



Gf->,r)=£ c 



hMr)| 2 

-e Q ±i7/2 ■ 



(11) 



In Eq. (|l 1[) , 7 is a level broadening. In an isolated sam- 
ple, one has to take the limit of 7 — > 0. Due to the 



discreteness of the spectrum of levels {e Q }, this extracts 
only e a the closest to the current energy slice e. Using 
the expression in Eq. (|lTJ) and taking the limit of 7 — > 0, 
one can so formalize the statistics of Eq. (^) that [gq 



^= A s(/^ ( ^ (r ' r ' ) - Gf(r ' r ' )) 



5(t-i<yVG? (r,r))> 



(12) 



The reformulation of Eq. ( p"2| ) in terms of the er-model 
exploits the fact [E2[ that any product of Green functions, 



1 [Gf(r,r)] n G^(r' ! r') = 



(13) 



sir' 



|s 2 (r)| 2 "e- L( * ) D*, 



can be represented as a functional integral over the 
8-component super-vector field * = (^1,^2), = 
^(Xm>Xm,C s m)' Tne super- vector * is composed 
of 4 anti-commuting and 4 commuting variables \ an d 
s, respectively. The indices m = 1,2 appear in order to 
distinguish between advanced and retarded Green func- 
tions. Besides "J, the charge-conjugate field ^> should 
be defined; one can find this definition in Ref. Q|. The 
action L(^f), 



£(¥)=»/ V(r)[e-H -U{r) 



i2A]#(r)dr, (14) 



A 11 = -A 22 = 1 

incorporates both the free-particle Hamiltonian H and 
the random impurity potential U(r). 

After Gaussian averaging over U |j, we derive a new 
Lagrangian with an interaction of the super-fields The 
interaction term can be decoupled by the functional in- 
tegration over a super-matrix field Q, so that any calcu- 
lation is finally reduced to the evaluation of a functional 
integral, 



/ DQ exp{-F[Q]}W(Q), 



(15) 



over the slow- varying superfields Q(r). This manipula- 
tion is analogous to the introduction of an effective order 
parameter in the theory of superconductivity. The free 
energy which determines weights of different configura- 
tions of Q appears after integrating over fast modes and 
has the form 

F[Q] = Jdr [-|Strln (-iH + f A + §) + f^StrQ 2 . 

(16) 

The 'anomalous mean' Q can be found from the self- 
consistency condition 
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>Q = J dp (-iH + 7A/2 + Q/2t 



(17) 



which minimizes .F[Q]. 

In the limit ol 7 — > 0, the solutions of Eq. ( |T7| ) take 
the values 



Q = VAV 



(18) 



from the degeneracy space of one of the graded symmetry 
group ]3q |. This field-theoretical model is strongly non- 
linear, since the matrix V satisfies the condition VV = 1 , 
and operations of the conjugation V — > V and the super- 
trace (Str) in Eq. ( |l~6|) are those defined in Ref. @ . The 
'rotation' V, 



1 - 1 u\ 





-iu 2 fu 2 







(19) 



is parametrized with the equal number of commuting 
and anti-commuting variables. The second matrix in the 
product in Eq. (|l^) is composed only of commuting ones. 
All anti-commuting variables are collected into the ma- 
trices Mi and v± . Smooth spatial variations of the Q-field 
at the length scale longer than the mean free path I and 
the influence of a small but finite value of 7 can be taken 
into account by the effective free energy functional, 



F[Q] = ^/drStr -f (VQ) 2 - 7AQ 



(20) 



The extension of this equation to the case of spin-i par- 
ticles can be found in Ref. §, too. 

To transform the formulae in Eq. (|[]|) into integrals 
over the Q-space, we expand the (5-function in Eq. Jl2] ) 
into the formal series in Gf~, and study the averages for 
all n 



[Gf(v,v)] n Jdr'Gt(v',v' 



Using Eq. Jl3|), each of them can be represented as a 
functional integral over the field ^ and then modified 
into the construction 



/..A-o J (2tt) 2 2(2n) 



In the latter equation, we perform an additional in- 
tegration over the phases £1,2 (which are hidden into 
the vectors v x = \/2(0, 0, e^ 1 , e"* Cl , 0, 0, 0, 0), v 2 = 
\Z2(0,0,0,0T),0,e 4C2 ,e- lC2 ) and add to the Lagrangian 
from Eq. (H) a weak perturbation 6L, 

5L= [ dr' [/i(ui*(r'))(*(r')«i) 



+A(« 2 *(r'))(*(r> 2 )5(r' - r)] 



After this, we have to evaluate the generating func- 
tional (/ D^e~ L ^- SL ^). Its exponent differs from 
that in Eqs. (fil]-E0h only by the perturbation 



SH = i / dr' [fivi <£> v± + \v 2 <8> V2&(y' — r)] 



added to the Hamiltonian Hq in Eq. (|16[). The latter 
results in an additional term 5F in the free energy func- 
tional; we find 5F by expanding the logarithmical expres- 
sion in Eq. ( |l6| ) into the series in [i and A. Doing that, 
we keep only the contributions up to the first order in /i, 
whereas 'cross-terms' which originate from pairing of 'f's 
at different coordinates (r and r') can be neglected. As 
an intermediate step, we obtain 

SF = ~ttv(viQvi) + ~ In [1 + \irv(v 2 Qv 2 )] , 

and after some algebra we arrive at 

f(t) =J§ lim^ / DQ [J * StrfaQtr))] (21) 



*(i-^Str(T(C)Q(r ))) e -*M 



where 



7T b 











and 



T(C)=7r 2 e^(ro+Ti)e-^. 



Everywhere below, t, are the Pauli matrices, and t is 
2x2 unit matrix. 



IV. REDUCED d-MODEL AND ITS 
SADDLE-POINT SOLUTIONS 

Basing on Eq. ( pij) , we can obtain the full statistics of 
local amplitudes for any regime. As we mentioned 
before, the universal expressions of Eqs. (||||) can be red- 
erived by assuming the zero-dimensional (0D) limit: the 
coordinate-dependent field Q(r) has to be replaced by its 
value at the observation point, Q(r ) = Q , which trans- 
forms the functional integral in Eq. ( ^l| ) into a definite 
integral over Q a . 

To go beyond the 0D approximation, one should take 
into account inhomogeneous fluctuations of the field Q. 
If we integrate over Q a = V AV , any of functions in Eqs. 
(HfJ) can be finally expressed in terms of relative rota- 
tions of the Q-field with respect to its value at r D . This is 
the reduced a-model. For its derivation, it is significant 
to note that the degeneracy space of the supermatrix Q 
is non-compact. Due to this property, the main contribu- 
tion to the integral in Eq. (|2l|) comes from the region of 



5 



large Q 's where Str(TQ D ) oc I/7. As a result, finite vari- 
ations of Q(r) produced by means of local rotations Q(r) 
— > Q(r') = V(r, r')Q(r)V(r, r') of the supermatrix field 
along the non-compact 'direction' can be taken into ac- 
count consistently, since they cover only relatively small 
environs of an 'infinitely large' Q D . 

Using the decomposition V (r) = V V (r), we define 
supermatrices Q of the reduced cr-model as 

Q = VAV, Q(r) = V Q(r)V , Q(r„) = A. (22) 

Due to the invariance of the Q~ s P ace i the transformation 
of Eq. ( p2| ) preserves the form of the gradient term in 
the free energy F in Eq. ((2^), whereas the second term 
in F can be modified as 



~- I <lr*U-[(l.Q[v)). 



Qo = V AV . 



A corresponding substitution can be done in the pre- 
exponential in Eq. j2l|), too. The explicit form and exact 
parametrization of the matrix Q a varies from one sym- 
metry ensemble to another. Nonetheless, in the limit of 
7^0, those parameters of the Q-space which are respon- 
sible for its non-compactness appear in the argument of 
the 5-function in Eq. ( pl|) in the same combination with 
the factor 7 that enters to the 'potential' part of the free 
energy, F2. Therefore, integrating over Q a in this limit, 
we eliminate 7 and convert Eq. ( |2l| ) into expressions 
which relate the distribution function f(t) to the gener- 
ating functionals represented in terms of the fields Q. 

The parametrization of Q-matrices, and, therefore, 
derivation and form of the reduced cr-model depend on 
the fundamental symmetry of the system. In parts A, B 
and C of this Section we specify this for unitary, orthog- 
onal and symplectic symmetry classes separately, though 
it turns out that the most essential part of our calcula- 
tion - the use of the saddle-point solution described in 
subsection IV. D - is quite similar for all of them. 



A. Unitary ensemble 

In the unitary case, the parametrization of the Q-field 
using Eq. (|l^ ) includes 'angles' 



(6t 

V 



< e < 7T 
< 6>i < 00 ' 



where only one of them is imaginary and makes the sym- 
metry group non-compact. Matrices U2 and V2 in Eq. 
( |l9| ) are diagonal and can be trivially eliminated from 
Eq. (|l]) as well as the external phases £. Other details 
of the integration over Q a are the same as those in Ref. 

H- 

The distribution function /, 



V dt 2 



(23) 



and the inverse participation numbers t n , n > 2, 

t n ee f °° t*f(t)dt = 2^ f °° t^<P u (t)dt, (24) 

can be related to the generating functional $ u (t) of the 
reduced cr-model, 

*»(*)= /g( ro)=A DQ(r) e -^[*.«]. (25) 

The free energy F u in Eq. ( |25| ) has the form 

F u [t, Q]=J drStr (VQ) 2 - \ AnQ(r)l , (26) 



and we remind that Q(r Q ) 
jection operator II in Eq. 



A at the origin. The pro- 
) is defined as 



7Tfc Tib 





TO 



(27) 



and selects from the Q-matrix only its non-compact sec- 
tor. 

The generating functional Q u (t) has several funny fea- 
tures. First, at t = 0, it has a completely invariant form, 
and, therefore, is equal to unity, what corresponds to the 
normalization of the wave functions, 



Vti = *„(0) = 1. 



(28) 



On the other hand, for any finite t, the reduced cr- 
model is a model with a broken symmetry, so that the 
free energy in Eq. ( p6| ) can be minimized by an inho- 
mogeneous solution Q(r). Indeed, |AII in the second 



term in Eq. ([26|) looks as a field tending to align the ma- 
trix Q along the non-compact 'direction' of the Q-space 
(related to the parameter 9\) 1 whereas the boundary con- 
dition at r D together with the gradient term is a rigidity 
attempting to prevent that. The competition between 
these two tendencies results in an optimal configuration 
of Q. To find such an optimal (saddle-point) solution, we 
use, again, the invariance of the Q-space with respect to 
rotations V and represent Q as 



Q(r) = Vt(r)Ai±gV t (r), P 



B 
B 



(29) 



where a weak perturbation P stands for fluctuations 
around the saddle-point, and the matrices B,B can be 
decomposed into blocks as follows, 



B = 



SXATQ + lSl t 2T3 

at 



S2,lT + is 2 ,2T 3 



(30) 



cr Q 
-at 



The form of the saddle-point Q = V t AVt, follows from 
the requirement of the absence of linear terms in the ex- 
pansion F u into the series 



G 



F u [t, Q]=F t + + F< 3 ) + F< 4 ) 
in the perturbation P. This selects 



\ ', = exp 











(31) 



(32) 



where the parameter 6 t (r) satisfies the optimum equation 



A0 t (r) 



txuD 



x 



(33) 



with the boundary conditions 6 t (r Q ) = in the origin and 
(nV)#t = at the surface of a sample. In Eq. ( ^3| ) A 
stands for the Laplacian in the real space. This equation 
is partly similar to the saddle-point equation derived by 
Muzykantskii and Khmclnitskii |17| when studying the 
problem of long-living current relaxation in open con- 
ductors, but it has a different non-linearity and - what is 
more important - different boundary conditions. 



B. Orthogonal ensemble 

The parametrization of Q-matrices in the orthogonal 
case is more complicated due to a larger number of inde- 
pendent parameters in it. In particular, the non-compact 
sector of the degeneracy space is parametrized by two 
'imaginary angles' - variables 6i t 2- 



< 6 < 7T, 
< 0i 2 < oo. 



0t o 
i{6iT + 6 2 n) 



Unitary matrices u 2 and v 2 in Eq. (|19|) have a more 
complicated form, too, 



M 
"■■2=1 e 1 ^ 3 J ' V2 



to 

e lXT3 



M 



1 — ifhf 
1 + ifhf 1 



where < 4>, \ < 27r, mi,2,3 are real, and the number of 
anti-commuting variables in m and V\ is twice as large 
as in the unitary case. 

The integration over Q can be performed in a com- 
plete analogy with the unitary case, but with several dis- 
tinguishing features. First of all, in the limit of 7 — > the 
main contribution comes from the region of the (61,62)- 
plane where cosh 6\ cosh 62 ~ 1/7- Since the product 
cosh 0i cosh 62 can be large at large 6\ as well as large 6 2 , 
we end up with the new form of projection operators, 



n. 



TTb(o) 
7Th(o) 



7Tb (o) 
7Tb (o) 



7Tb (o) 





TQ + Tl 



in the free energy 



F [t', Q] = J drStr ^(VQ) 2 - ^An o Q(r)j . (34) 

This n o determines the direction of an effective 'force' 
along the symmetrically chosen non-compact 'direction' 
(61+62). 



Next, in the orthogonal ensemble, one has to keep the 
external phase factor e'^ in Eq. (gl|) until the end of 
the integration over Q a , which results in the integro- 
diffcrcntial relation 



fo{t) = virS{^ d ^o(t + ^)} 



(35) 



between the distribution function f Q and generating func- 
tional 



<W = /Q (ro)=A DQexp(-F [<',Q]). 



(36) 



The generating functional ^ (t') gives directly the inverse 
participation numbers t n , n > 2, 

tn = ^ TJ ^^(tT- 2 *o(t')dt>, (37) 

and, again, has the property $ o (0) = Vt\ = 1. 

To study the fluctuations near the saddle-point, we 
represent Q in the form of Eq. (E9f) , where 



B = 



B- 



V2 



B+ = 



_ ( (s u t + is^ 2 T 3 )r± 



(4 ±} ) 



(S^T + isf 2 T 3 )T± 



and 



r (±) 



±\* ) , T± = T ± Tl. 



In this decomposition, s a g are real numbers, and 
(&&)* ~ anti-commuting variables. Indices '±' are intro- 
duced for the later convenience. Everywhere below, we 
keep superscript '-' but omit '+'. 

Similarly to the unitary case, the free energy F a which 
governs the generating functional $ G has the minimum 
at Q = V t AV t , 



V t = exp 







± 4 




where the variables 6 t = 6i + 6 2 , 6 t = 6\ — 6 2 and x, cf> 
satisfy the following equations: 

A6 t (v) = -^e- 9 *, Afl t "(r) = 0, e*^ = -1. (38) 

The conditions at the origin and boundary are the same, 
as in the unitary case. The latter gives 6j (r) = 0, and the 
non-trivial saddle-point is related only to the symmetric 
variable 8 t . 
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C. Symplectic ensemble 



and 



An analogous investigation of the statistics of spin- 
polarized electron waves in the case of a strong spin-flip 
scattering needs an extension of dimensions of Q-matrices 
and the following analysis of the degeneracy space related 
to their gapless Goldstone modes. The gaps in the spec- 
trum of Q's appear due to a large spin-relaxation rate, 
t.7 1 , which can be caused both by the spin-orbit coupling 
built into the material properties or by the spin-flip scat- 
tering on a classical randomly oriented static magnetic 
field. In the former case, the time-reversal symmetry is 
conserved, while in the latter this invariance is violated 
by the source of a scattering. Since triplet components of 
Q correspond to gapfull modes [§,0, only singlet modes 
have to be taken into account, so that the number of in- 
dependent variables of Q is the same as in the spinless 
case. One has to remember only that all the matrix ele- 
ments of Q are multiplied by the unit 2x2 spin-matrix 

To- 
lll this subsection, we work with the distribution /(J., t) 
of a local spin density t = \ip± | 2 of the eigenstates defined 
by Eg. (Q) and above it. An incorporation of spins into 
Eq. ( |2l| ) can be done by substituting 

7Tl,2 — > 7Tl,2 ® To, T — > T ® f|, 

where stands for the direct product of matrices, and 

1 0\ . 
x 1 J' Ti = 

0' 
1 

along a single direction, and the integration over Q Q gives 
us 

/a><)=2F%^ $ s (t)=/DQ(r)e-^[*.«], (39) 
where 



?i are the Pauli spin operators: fo = 

The degeneracy space of Q Q is non-compact 



F s [t, Q]=J drStr ^f- (VQ) 2 - \ AITQ(r) 



(40) 



and II is exactly the same as in Eq. (p7|). 

From the point of view of the rest of calculations, the 
case with a broken time-reversal invariance is equivalent 
to the spinless unitary symmetry class 0. That's why 
we mark the quantities related to this symmetry with a 
label 'sit' and generate the distribution function f su {l,t) 
form the distribution function f u (t) in Eqs. j23]-p6|) as 

f su (i,t) = 2f u (2t)atD^2D. 

As concerns the symplectic ensemble with the time- 
reversal symmetry, it demands an extra calculation, since 
in the presence of the spin-orbit interactions it cannot be 
reduced to the spinless orthogonal one. The parametriza- 
tion of the Q-space in this case is given by Eqs. ( pj| ) with 



/ 8t + 6'n 
^ i6 lTo 



< 9,9' < 7T, 
< 6»i < oo, 



"2 



M 



, V2 



e tXT3 \ M _ 1 - irfiT 
r / ' 1 + im? 



where < <f>, \ < 2n, and mi. 2, 3 are real variables. 
After this, the saddle-point configuration of Q for the 
symplectic-orthogonal case can be found as Q = V t AV t , 



V, = exp 







'fM 




where 9t satisfies the optimum equation and the matrix 
M is chosen in such a way that M = — tq: 



A0 t (r) 



00, 



(41) 



with the boundary conditions 9 t (v ) = and (riS/)9 t = 
at the surface. 



D. Optimal free energy and fluctuations near the 
saddle-point 

After comparing the saddle-point equations in Eqs. 
([33|,[58],[ll]) , one finds that they are similar in different 
symmetry classes. The difference between the unitary, 
orthogonal and symplectic ensembles leads only to dif- 
ferent values of a coefficient /?, 



(42) 



in the expression for the optimal free energy 

F t =f3fdr{^(V9 t ) 2 +te- e t}, (43) 

and in the higher order terms of the expansion of F[t, Q] 
in the environs of the saddlcvpoint. The generating func- 
tional $(t) from Eqs. ( f2l|[36| , [39"| ) can be represented in 
the form 

$(*) = J(*)exp(-F t ). 

In the conducting regime, the value of the optimal free 
energy determines the leading term in the exponential of 
the generating functional, whereas the effect of fluctua- 
tions around the saddle-point is included into the func- 
tion J(t), 



J(t) =/DPexp{-F( 2 ) -F<< 4 \..}. 



(44) 



Due to the normalization condition in Eq. (|28|) , the re- 
lation J(0) = 1 holds exactly, and the contribution from 
the fluctuations P can be calculated by expanding the ex- 
ponential in the integrand in Eq. ( pi| ) into the series in 
the higher-order terms F ( 3 > 4 >---) and performing Gaussian 
integrations over with the weight exp{— F^} determined 
by the second-order correction to the free energy. 

The applicability of such a perturbation theory is justi- 
fied by the fact that the higher orders are, at least, by the 
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factor of (2ttvD) 1 <C 1 smaller, as compared to what is 
given by 



J(t) » /DPexp{-F( 2 )[<,P]}. 



(45) 



The latter is nothing but the super-determinant of 
the Hamiltonian related to the fluctuations around the 
saddle-point. The value of J(t) differs from unity merely 
because the symmetry between fermionic and bosonic de- 
grees of freedom is broken by the optimal solution. Since 
not all the projections of the infinitesimal P to the gen- 
erators of the Lie algebra of the graded symmetry group 
arc cquivalcntly affected by the symmetry breaking, it is 
convenient to separate in the terms which feel the 
existence of the optimal solution from those which do 
not. Depending on the physical symmetry class, this in- 
volves different sets of variables. Nevertheless, after an 
appropriate diagonalization, quadratic form can be 
represented uniquely for all symmetry classes: 



F, 



(2) 



(2) 



Su,SU 







(46) 



The term Fq 2 ^ in Eq. (|46|) is composed of fields which 
are not affected by the symmetry breaking, 

F (2) = 2-kvD Wdr [8az8{a~y + &-£(*-)*}■ 

a=l,2 

This term does not contribute to the function J in the 
Gaussian approximation, due to the symmetry between 
boson and fermion degrees of freedom incorporated in 
it. On the contrary, the first term in Eq.([i6|) is the sum 
over those four pairs of dynamical variables which feel 
the violation of the boson-fermion symmetry 

F t {2) = 2-kvD I dr i J2 fivcxdal + U a a a a* a ] 

la=l,2 



Due to that, the pre-exponcntial J can be represented as 
J = exp { | E„ In ((xAn))V Yl a<f} X?») } , (47) 

where the sum is extended over all the eigenvalues of the 
spectral problem 

[-A + U- X ]<f> = 0, 0(r o ) = 0, nV^(S) = 0. (48) 

As we already mentioned above, in the quadratic ap- 
proximation, any difference of J from unity is due to the 
broken symmetry between fermionic and bosonic degrees 
of freedom in P. The broken symmetry in the Hamil- 

(2) 

tonian F t is the consequence of the difference between 
the effective potentials 



(49) 



in the Hamiltonian F^ 2 \ In Eq. (|49|), k a \ = g Q i = 
and ki.2 = 4, fc 2 ,2 = 0, q a ,2 = 2. The spectra {Xa(n)} 
and {Xs( n )} °f modes remain gapfull, since all U > 0. 
Moreover, due to the sum rule 



(50) 



their main contribution to J comes from low-lying eigen- 
values of Eq. (flsl) . Since the set of IPs in Eq. (|49| ) 
depends on the optimal solution alone, the calculation of 
the correction to the exponent in this order in (2ttvD)~ 1 
can be performed simultaneously for all symmetry classes 
and is small. 

On the other hand, the effect of fluctuations can be- 
come important once we want to extend the considera- 
tion of samples with the size larger than the localization 
length. Such a step which is left beyond of the scope 
of this paper would need some kind of a renormalization 
of the saddle-point. We would only like to stress that 
this can be a way to avoid the previously found difficulty 
p6|-|38| to obtain a stable fixed-point in the renormaliza- 
tion group treatment of the localization problem. 

The existence of the saddle-point and relatively small 
contribution from fluctuations in the metallic regime 
makes it easy to find the form of the cross-correlation 
function R(t,r) from Eq. (||), too. If we study the enve- 
lope of the wave function at large enough distances r ^> I 
from the position of a high amplitude splash, the reason- 
ing used above can be repeated for R(t, r) with minor 
modifications, and we arrive at 



R(t, r)/f(t) cx te- 0tir \ 



(51) 



so that one can say that the envelope of |V , ( r )| 2 follows 
the form of the saddle-point configuration of the reduced 
cr-model. 

Up to now, we performed our analysis without refer- 
ring to any specific dimensionality of the system. On the 
other hand, from Eqs. d3^f^ , f4f| , [48| ) , one can see that 
the saddle-point solution, and, therefore, the optimal free 
energy F t crucially depend on the dimensionality. In the 
next Sections, we discuss the statistics of local ampli- 
tudes of wave functions in ID, 3D and 2D conductors 
separately. 



V. NEARLY LOCALIZED STATES IN A SHORT 
DISORDERED WIRE 

It is well known that the localization effects are the 
strongest in one- and quasi-one-dimensional (ID) con- 
ductors |^|,^). Even when disorder is weak, the quan- 
tum diffusion of a particle is blocked at the length scale 
L c = [32iri/D, where D is the classical diffusion coefficient 
determined by the elastic impurity scattering M and (3 
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is specified by Eq. (|42|). Since in the quasi-lD wire the 
effective density of states v is that integrated over cross- 
sectional width w or an area, w 2 , the localization length 
L c ~ l(w I '\F) d ~ x can be much longer than the mean free 
path I. This allows us to consider the short wires L < L c 
with an already developed diffusive regime, and below we 
analyze the distribution of amplitudes and the shape of 
untypical states which are the precursors of localization 
at larger distances. Since there is a lot of known about 
the ID case [p p9|j40| , its example can be a good point 
to compare the results of the saddle-point approach with 
calculations based on the use of the exact transfer-matrix 
method Q|. In fact, in the ID case, the saddle-point cal- 
culation is nothing but the 'semiclassicaP solution of the 
effective Schroedinger equation on the Q-space JtJ which 
appears in the transfer-matrix method. Such a 'semiclas- 
sics' not only gives the results which are very close to the 
exact solutions (io| , but also enables us to make a state- 
ment about the form of typical pre-localized states in the 
metallic regime. 

In the following paragraphs, we apply the scheme of 
calculus described in the previous section. First of all, 
we have to solve the saddle-point equation, 



d 2 Mx) 



t 

■kuD 



exp{-0 f }, 



(52) 



and use its solution 9t(x) when calculating the optimal 
free energy F t . Due to the condition 6t(x ) = at the 
observation point x a , the latter splits the wire with the 
length L into two (not necessarily equal) intervals < 
x < Li=L,R- The form of 9t{x) can be found separately in 
each of them. In dimension one, the differential equation 
in Eq. ( |52] ) can be solved exactly and we represent 
its general solution in the form 



x > 0; i = L, R. (53) 



s{A lV /T7(l-^)} 



of 



Although one can notice that the general form 
this solution formally contains a singularity at = 
—Li[ir I '(2Ai\fTi) — 1], the latter is illusory since it takes 
place in the non-physical region x < and plays no role 
unless it comes up to the formulation of the limitations 
to our theory. The requirement Xoo 3> I which emerges 
from the existence of the singular point is related to the 
conditions on maximal values of gradients permitted by 
the use of only the lowest-order gradient expansion terms 
in the free energy functional in Eq. (|l]). We shall dis- 
cuss the consequences of this condition at the end of the 
Section, assuming for a while that it is satisfied. If so, the 
consistency equations on the parameters Ai come from 
the boundary conditions d x 8t(Li) — at the edges and 
have the form 



Ai = cos (A l ^/T~) 
where T are defined as 



Ti = 



2nvD ■ 



(54) 



(55) 



The optimal value of the free energy can be found, in its 
turn, as 



(56) 



In general, the exact form of F t in Eq. ( |56| ) based on 
the closed set of equations in Eqs. ( |53U5q ) can be studied 
numericaly at any values of the parameters included, but 
a somewhat simpler analytical expression can be written 
in the asymptotical regions. First of all, we examine 
the limit of small amplitudes, Ti < 1, where the exact 
distribution has to match with the random matrix theory 
results. At T t < 1, the results of Eqs. (|5|-|6|) can be 
expanded into the series in Ti, which gives Ai w 1 — 
and 



IT 

2 1 



llT 2 

24 i 



F t w Vt 1 - 



TU 



^ 3L 

i=L,R 



We see that in this limit the Porter-Thomas formulae, 
Eq. (^8]) give a good approximation for the body of the 
distribution function f(t) which describes those ampli- 
tudes t which are t < L~ l ^/L c /L. Otherwise, the second 
term of this expansion, VtTi, becomes larger than unity 
and strongly affects the probability to detect a too high 
splash of the wave function. 

When Ti S> 1, the optimal c onfigura tion e~ 0t develops 
at the length scale of A = w2TwD/t, where it can be 
approximated as 



-9 t (x) 



{\/xf 



(57) 



and gets less sensitive to the presence of boundaries. In- 
deed, in the limit of T ^> 1, one has Ai 



-1/2 



§77 1/2 (i 



.), and the exact expression for the optimal free 



energy can be expanded into the series in T 



-1/2 



Sr}- 



(58) 



where 



\T 



-1/2 



T" 

2 1 



i = L,R. 



The leading term in Eq. ( |5q ) does not depend on the 
system length and (in the orthogonal ensemble) coincides 
with the asymptotical behavior of the distribution func- 
tion of local amplitudes in infinite wires. The latter has 
been found by Mirlin and Fyodorov on the basis of 
the analysis of the transfer-matrix equations derived by 
Efetov and Larkin |39| . Although we did our calculations 
for the short-length samples, L < L c = j32nvD, the re- 
sults for the high amplitude splashes surprisingly agree 
with those for the infinite geometry, even up to the lead- 
ing term of pre-exponential factor J. The latter is derived 
in Appendix A by taking into account fluctuations near 
the saddle-point. Roughly speaking, the 'semiclassical' 
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solution of the transfer-matrix equation gives an almost 
exact result. The full form of the tails of f(t) at t > g/L, 
where g = L c /L, can be represented as 



/W(t)«C 



exp [-A*J]JLJ{1 



(59) 



and is applicable up to the amplitudes t ~ L c /l 2 ~ 
w/(1\f)- The latter condition emerges from the require- 
ment of a smoothness of the saddle-point solution, so that 
its characteristic length scale A should be longer than the 
mean free path, A > I. Otherwise, the singularity of the 
equation in Eq. (|53) comes too close to the physical 
part of the space (x > 0) which would create too large 
gradients forbidden within the framework of the used ap- 
proximations. 

The distribution function given by Eq. (|59| ) evidences 
that the states which are responsible for the rare event we 
discussed in the previous paragraph are (at least, partly) 
localized. Nevertheless, even for the largest amplitudes 
t > L C J L 2 , the effect of edges is still present, which means 
that this in not an exponential localization. On the ba- 
sis of an analysis of the cross-correlation function R(t, r) 
from Eqs. (§) and @, we can say that, within the range 
of distances x < L c , the envelope of the density of pre- 
localized states ipt(x) resembles the form of the optimal 
solution and has the power-law intermediate asymptotics 



(60) 



In contrast to the 2D case which we discuss in Section 
VII, the derived form of a typical wave function has the 
same exponent for all pre-localized states, independently 
of the amplitude of their top splashes. Further, the tails 
of the envelope in Eq. (J60Q are integrable, so that the 
inverse participation numbers which one can find on the 
basis of Eqs. (|) and ( pi] , [37| ) do not indicate any fractal 
behavior. 



VI. EIGENSTATES STATISTICS IN D = 3 

The localization effects in weakly disordered 3D con- 
ductors are known to be the least pronounced [9 , so that 
the eigenstates statistics of 3D conductors has to be most 
similar to the universal one. Nevertheless, even in this 
case, not all states are described by the Porter-Thomas 
distribution, and this Section is devoted to the disorder- 
dependent corrections to formulae in Eqs. in di- 
mension three. 

To describe the statistics beyond the universality limit, 
we have to solve optimum equations (|33| , |3"8| , [4l|) . For the 
sake of simplicity, we consider a spherically shaped con- 
ducting particle and place the observation point into its 
center. This enables us to seek for the solutions 9 t (r) in 
a symmetric form. Nonetheless, even that does not help 
us to find the exact form j4^] of a general solution of the 
non-linear equation 



A9 t (r) = [r- 2 d r r 2 d r ]e t 



(61) 



so that we have to develop the following approximate 
procedure. The non-linear cr-model, Eq. (p0|), was de- 
rived under the conjecture of smoothly varying Q-fields. 
This implies that the distances shorter than the mean free 
path I are excluded from our consideration, and the con- 
dition 9t(0) — at the origin has to be substituted by the 
condition 9 t (ro) = at the sphere of a radius Tq ~ I. Af- 
ter this, we scale the distances by I, so that 9 t = 9 t (r/l), 
and solve the problem iteratively. The iterative proce- 
dure appeals to the fact that the Laplace equation which 
one can get by neglecting the right hand side of Eq. ( |6l| ) 
has non-zero solutions and that the parameter 



P = ^TTvD/(l 2 t) 

which appears after rescaling the distances r with the 
mean free path is large. The latter condition restricts 
our consideration to the amplitudes t < (ZA^) -1 which 
are smaller than the density formed by the forward-and- 
backward scattered trajectory between two impurities. 
As the first step, we expand the desired function 9t{r) 

as 



l(0) 



9[ 0) =A(l-l/r). 



The first term in it satisfies the Laplace equation but does 
not satisfy the necessary boundary condition at r = L. 
The term 9^ is added in order to satisfy the requirement 
d r 9 t = at the external edge. It must turn to zero at 
r = I and can be found from the linearized equation 



\u- 2 d u u 2 d u ] 



! exp(-0 t (O) ), 



9 t (1) (l) = 0, 



where u — r/l. After this, the non-linearity of Eq. (|6l| ) 
transforms into a self-consistent determination of the pa- 
rameter A from the algebraic equation 



A = P - 



L/l 



1 



2 dwexp{-A(l-u- 1 )} pb — - 



3p 2 



L 



which arises from the requirement d u 6^ (L/l) = 0. One 
could continue the iterative scheme even further and add 
corrections which improve the function 9^ (x) itself, and 
so on, but this is not necessary for evaluating the leading 
terms of the optimal free energy F t . So, we stop the 
iteration after the first step and find that 



exp(-A(l~l/r)), r>l. 



(62) 



The combination of the parameters which stands in 
the right-hand side of Eq. (^lj) and the self-consistency 
equation itself can be rewritten in the form 



Ae A = T = ^% 



Vt/(p F l) 2 



(63) 



where 773 ~ (2/) 1 and the condition ppl 3> 1 corre- 
sponds to the limit of a weak disorder. The optimal free 
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energy related to the saddle-point configuration can be 
calculated, too, and has the form 



+ 4} 



(64) 



When T <C 1, the calculation both of the self- 
consistent A and the related value of F t can be performed 
as an expansion into a series in the parameter T, i.e., we 
approximate A rs T — T 2 + ... and 



F t = (3Vt 



1 - 



±ll 2 vD 



m ~ (20" 



When T ^> 1, the leading terms arise from the es- 
timation A ~ InT. In each of these two regimes, the 
generating functional <5(i) has the form 




ln 2 T- 



, T< 1, 
}, T>1, 



(65) 



which can be used for evaluating both the distribu- 
tion function /(*) and IPN's using Eqs. and 
(|J,^). At this point, we have to mention that the coef- 
ficient r]3 ~ (2Z) _1 cannot be determined better than by 
the order of magnitude. We also remind that different 
symmetry classes are distinguished by the parameter (3: 

Eq. ( pq ) indicates that the noticeable deviations from 
the universal Porter-Thomas distribution start rising at 
local densities tV ~ PfI (the second term in the expo- 
nent in Eq. ( |65| ) becomes larger than unity) and then 
develop into the logarithmically-normal asymptotics at 
tV > (pfI) 2 - On one hand the states which generate 
such an asymptotical tail are not typically metallic. On 
the other hand, both the form of the envelope, 



|V>t(r)| 2 oc exp{-A(l - i/r)}, 



(66) 



which we extract from the shape of the optimal solution 
in Eq. ( |62"| ) and the scaling of IPN's with the integer 
power of the system volume for any n, 



min{y>(n), l2ir 2 isD/r] 3 ] n } f n 2 rj 3 



yn 



/34ir 2 vD 



<p u = n\, tp a = (2n - 1)!!, y s = n\/2 n , 

indicate that these 3D states are not localized in a con- 
ventional sense: They always have a finite part of the 
density 'equally' smeared all over the sample. Of course, 
these extended density-tails decrease when t approaches 
the limiting value t ~ (IXp)^ 1 , but our methods do not 
allow us to make a statement about the structure of 
standing waves at the scale of r < I. 

The version of the a- model we used above also restricts 
our consideration to the metallic regime pf I S> 1. The 
development of a theory at critical conditions, ppl < 1, 



requires the use of more sophisticated methods p|[ . Nev- 
ertheless, the common believe which arises from the most 
of known localization theories || is that the dimension 
d = 2 is critical, so that the analysis of wave function 
statistics in 2D disordered conductors would manifest the 
important features of the criticality. 



VII. MULTIFRACTALITY OF EIGENSTATES IN 
WEAKLY DISORDERED 2D CONDUCTORS 

To find an optimal configuration in the 2D case we 
limit the length scale of the Q-fidd variations from be- 
low by the value of a mean free path - similarly to what 
we discussed in the previous Section. This modifies the 
boundary conditions into Otifo) = at vq ~ /. For the 
sake of simplicity, we consider the sample in the form 
of a disk (with the radius L) and place the observation 
point r„ right in its center. Then, we seek for an axially 
symmetric solution 0t (r) of the equation 



(r^drrd^Ot = - 



t 



■kvD 



A. Exact solution 

This can be done both exactly [Q or using an iterative 
procedure developed in Section VI. The exact solution 
can be represented in the form 



-9 t _ 



(67) 



where p 



2 -jfp-, and A has to be found from the 
boundary condition at the sample edge r — L, 



^ T T P ^ + A = ^J^ 



(68) 



After substituting the saddle-point solution 9 t (r) from 
Eqs. ( |67| to Eq. (^3|), we also find the optimal free 
energy, 



F t = (3An 2 vD I In 



( (L/l) 



(l+A 2 ) 



\ p 2 [l - A 2 ] 



+ 2{l~y/A 2 +p- 2 ) 



Together with Eq. (p8|) , the latter expression can be 
used for the numerical analysis of the exponential of the 
distribution function. 

The numerical analysis shows that the consistency 
equation in Eq. ( |68|) has positive roots only if p > 
ln(L/Z) 3> 1, which provides a reasonable limitation to 
the wave functions amplitudes we can study using this 
method: We have to restrict the density t of a splash by 
the value (A^ /) _1 related to the density of an electron 
state bound to the forward-and-backward scattered tra- 
jectory between two impurities. At the same time, in the 
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limit of p ^> 1, the roots of Eq. ( p8[ ) can be approximated 
by A = 1 — fi, where fi < 1. The same conditions gives 
us a possibility to replace the exact solution in Eq. J67| ) 
by 



(l/r) 



2\i 



what means that there is an easier way to get a satisfac- 
tory approximate solution of the saddle-point equation 
in d = 2 similar to what we did in d = 3. 



B. Solution using iterations 

The result of Eq. ( |69| ) can also be derived using the 
iterative scheme discussed in the Section VI. Being ap- 
proved by the strong inequality /) » 1, we solve, first, 
the linear Laplace equation by choosing its solution in 
the form which satisfies the boundary conditions at the 
origin, 



3 (0) 



= 2^1n(r/0, 



where the parameter fj, will be the subject of the next 
iteration. That's, we seek for such 9 t = 9^ + 9^(r/l) 
which satisfies the boundary condition at the external 
edge, and where 



[u- 1 d u ud u }9[ 1) 



-2p~ 2 u 



-2p, 



9 t (1) (l) = 0, u = r/l. 



This gives 0^' in the form 



3(1)/ 



2(1 -fif 



1 



u 



2(1-M) 



The requirement d r 9 t (L) — gives rise to the consistency 
equation and enables us to formulate the approximate 
procedure in the closed form, 



2/J, : 



(L/l) 



2-2a 



(1 - ■ 



The use of the iterative procedure is formally limited by 
a requirement 0[ <C 1. 

The parameter [i in the above equation can be found 
(with the accuracy controlled by l/ln(L/Z) <C 1) as 



21n(L/i) ' 



ze 



T 



tV\n(L/l) 
2ix' i vD 



(70) 



and varies when the amplitude t changes. For example, 
the cross-over of the optimal solution to the homogeneous 
Q = A occurs in the limit of T <C 1 where one can find 
that 

\t{1 - T +...)/ In j 
In the opposite limit of large amplitudes, T 1, 
/i«llnT/lnj<l. 



The approximate form of the optimal free energy can 
be found, in its turn, as 



F t fa /34tt 2 vD{[i + p 2 In f }. 



(71) 



(69) When T <C 1, this can be expanded in T as 



F t « pvt{l - | + ...}. 

When T > 1 (but still t < (/A F ) _1 ), the leading term in 
the optimal Ft takes the form 



F t ps P-k 2 vD 



InT 
\n{L/l)' 



Although the size of the system enters these formulae, 
the logarithmically weak dependence of F t on L makes it 
meaningful to use the derived expressions for an arbitrary 
position of the observation point inside the sample of an 
arbitrary convex shape. 

The fluctuations around the saddle-point configuration 
and the resulting pre-exponential factor J(t) for the 2D 
case are discussed in the Appendix B. All over the con- 
duction regime, their contribution to the value of the 
generating functional is small, as compared to that of 
the optimal solution itself. 



C. Distribution function and IPN's 



All this enables us to calculate the distribution func- 
tion /. For small amplitudes t < 2-kvD/[L 2 In j], one 
obtains 



«Fexp {-pVt [1-f 




1, U > ( 72 ) 

s; su 



where T is defined in Eq. (|7^), and 

Po — 2' @ u — 1' Ps.su 2. 

In the opposite limit, t > 2nvD/{L 2 In j], the distribu- 
tion function takes the form 



/ (2) ^ex P (-g-ln 2 T 



(73) 



Eqs. ((72]) and ( |73| ) generalize our earlier result ^] to var- 
ious symmetry classes. They show that for any of the fun- 
damental symmetries - orthogonal, symplectic and uni- 
tary - disorder makes the appearance of high-amplitude 
splashes of wave functions much more probable than one 
would expect from the Porter-Thomas formula and, as 
concerns the most extra-ordinary events, tends the tails 
to take the logarithmically-normal form. When being 
written for the orthogonal ensemble, the logarithmically- 
normal law in Eq. (pq) strikingly coincides with the form 
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of the asymptotics of the distribution of the local density 
of states and conductance fluctuations in open systems 
found in Ref . || , although our theory was made for closed 
systems and is based on a different scheme of calculations. 
This agreement reveals the deep relationship of these two 
results obviously caused by the localization effects. But 
the localization of wave functions which are responsible 
for the asymptotic events is not the localization of a par- 
ticle in a confining potential. The tails of these states 
do not decay exponentially: Even in the asymptotical 
regime T > 1, the size L of the system influences the 
distribution. The splashes look as if they were formed by 
focusing the waves by some rare configurations of scat- 
terers. The structure of these states can be anticipated 
from the way how their distribution feels the boundary 
or - directl y - from the cross-correlation functions R(t, r) 
in Eqs. (| 51). Following the form of the optimal solu- 
tion, the envelope of the density of such a state has a 
power-law asymptotic tail 



\M*)\' 



(l/r) 



2/i 



(74) 



which approaches the limiting r~ 2 dependence for the 
highest amplitudes t ~ (ZAj?) -1 . 

Moreover, the form of IPN's, t n derived on the basis 
of Eqs. fl2^j37| ) shows such a scaling with the size of a 
system which implies them a multifractal structure. To 
find the moments t n accurately enough, we have to take 
into account that, although the cross-over to the OD case 
looks like the formal limit T(t) — > 0, the Porter-Thomas 
statistics fails unless the condition tV <C \T2mvD is sat- 
isfied (see Eq. (frgl)). Hence, only first few ratios t n , 



2 < n <C v2vn7D, can be estimated using a finite poly- 
nomial expansion of f(t) into the series in T, and their 
first terms reproduce corrections to the universal statis- 
tics found pertubatively in Ref. [^7|. Alternatively, we 
derive the higher order IPN's from Eqs. fl |2~i| , |36| ) using 
the saddle-point method. The moments t n calculated in 
both ways are in a good agreement with each other and, 
in the leading order, can be represented as 



t, 



where 



min{ip(7i) ,[2ttuD / In i]"} 

Wv 



/3~ 



(*) 



(n-l)d*/2 



4tt 2 vD ' 



(75) 



(76) 



As one can see from Eqs. ( |75| , |76| ) , we end up with 
such a volume-dependence of the inverse participation 
numbers t n that manifests the multifractal behavior of 
quantum states, Eq. (^|). The multifractality seems to 
be the generic property of 2D disordered systems. The 
multifractal dimensions in Eq. (76) are calculated in the 
leading order on the inverse conductivity, so that all over 
the metallic regime, the dependence of d* on n and disor- 
der is accurate enough and qualitatively agrees with nu- 
merical results ]ll]] . Due to the limitation t < (IXp)^ 1 , 
the above equations work at n < 2nvD, so that n— 5 > 0, 
and the fractal dimensions d* in Eq. (f76|) are positive. 



VIII. DISCUSSIONS 

Summarizing the results of the paper, we studied the 
manifestation of precursors of localization among the 
eigenstates of isolated disordered conductors with the size 
smaller than the localization length. In order to detect 
them, we analyzed the statistics of local amplitudes of 
wave functions, t = |-0| 2 , and, at the tails t ^> V^ 1 found 
strong deviations from the universal Porter-Thomas dis- 
tribution (see Eqs. (||-||)) associated with the typically 
extended-type behavior. The universal statistics equally 
describes the quantum states of various classically chaotic 
systems; it depends on their fundamental symmetry but 
not on the physical dimensionality or the level of disor- 
der. Such a description can be successfully applied to 
the most of the states (extended ones) in the metallic 
regime and gives the body of the distribution function 
of their local amplitudes. The deviations from the uni- 
versal laws start rising at the amplitudes t ~ \fg/V and 
finally develop into a completely different asymptotics at 
t ~ g/V. In dimensions d = 2 and 3, the form of the 
asymptotics is described by the logarithmically-normal 
tails in Eq. ( |To| , |65| , |73| ) , whereas in d = 1 it has a stronger 
dependence, / cx exp{— Ay//3tL c }. 

The scheme of calculus (see Sections III and IV) and 
the similarity between our results for isolated systems 
and the asymptotics of distributions of the local density 
of states and fluctuations of other quantities in ID pQ ] 
and 2D conductors indicate that the above-mentioned 
long tails are strongly influenced by the localization. 
To answer the question, how the localization develops, 
we can refer to the fact that the deviations from the 
Porter-Thomas distribution appear as a so small number 
of events, oc exp(— y/g), that their occurrence near the 
Fermi level in a specific sample is a typically mesoscopic 
phenomenon. We interpret this as that the rare top- 
amplitude splashes are not locally implicit as portions to 
any state but represent very non-trivial configurations of 
waves which can be found more and more often if the 
disorder increases. The analysis of the cross-correlations 
R(t, r) also indicates that the states which are responsible 
for locally the highest amplitudes ^(fo)! 2 > g/V have 
individually specific envelopes of their decaying density 
far away from the observation point r Q . In ID and 2D 
systems, the tails of the envelopes obey the power-law 
dependence \tp(r + r)| 2 oc r~ 2fl . In d = 1, its expo- 
nent is /j, — 1, so that the density of these tails is per- 
fectly integrable at long distances, and one could speak 
about them as about nearly localized ones. In dimen- 
sion d = 2 the exponent fj, is limited by fj,(t) < 1 (so 
that it is not the localization in the usual sense) and is 
individual for the states with different amplitudes of the 
top-amplitude splash. Such a behavior of pre-localized 
states in d — 2 coexists with the typically multifractal 
behavior of the inverse participation numbers which has 
been observed earlier in various numerical simulations at 
the critical conditions of the localization-delocalization 
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transition |1Q-|15|]. Unfortunately, at the present stage 
we cannot approach close to the 3D Anderson transition 
and to check the multifractality globally. Nevertheless, 
even in d = 3 we find the non-trivial logarithmically- 
normal asymptotic behavior of the statistics, although 
the states which seem to be responsible for that are not 
localized. 

The combination of the facts presented above forces us 
to suggest that the details of the structure and unusual 
statistics of rare pre-localized states which we discussed 
in the present paper have something to do with the statis- 
tics of extraordinary multiply self-crossing diffusive paths 
which anomalously often return to the same spatial co- 
ordinates r Q . This means that the almost confined states 
develop because of rare shortened classical trajectories 
which form a closed loop not only in the real space, but 
also in the full phase space, since they finally come to 
the same 'sell' dpdx = Ti where they started. For in- 
stance, the most dense configuration could be formed by 
a forward-and-backward scattered trajectory of a particle 
bouncing few impurities. An additional argument which 
supports this scenario relates to the conditions limiting 
the validity of our theory. Basing on the use of the a- 
model, we have necessarily to cut the linear length scale 
of the supermatrix Q-field variation from below by the 
value of the mean free path. Nevertheless, the densities 
which can be described in our approach are limited by 
the value l/(7Ap -1 ) instead of a naively expected inverse 
volume l~ d . This is only possible if the states which we 
study are locally anysotropic at the fine scale of distances 
of about I and typically have a snake-like structure with 
the cross-sectional width ~ Xir 1 . 



with lengths Ll,r, where each of Jl,r is determine d by 
the eigenvalues of the Schroedinger equation in Eq. ( ]4q ) 
with the symmetry-breaking potentials 



U, 



kT, 



I sin 



1 



where i = L,R. When T <C 1 (in the paragraph below, 
we omit indices L and R), these potentials can be treated 
perturbatively. Their first-order corrections cancel due 
to the sum rule from Eq. (po|), so that J » 1 + T 2 w 1. 
When T> 1, the same cancellation eliminates contribu- 
tions from the high-excitation eigenvalues \ > ( 7r /2i) 2 7 1 , 
so that the important contribution comes from the low- 
energy part of the spectrum, x < /2L) 2 T ', where one 
can approximate 

Urn (7r/2L) 2 [k + k'/sin 2 (nx/2L)}, k' = k + n. 

Using this approximation, the spectral problem of ID Eq. 
(E^) can be solved exactly. To find the exact solution, 
one has to change variables from x to y = cot(7rx/2L) 
and then seek for solutions in the form <p = P n {y 2 )/(1 + 
y 2 ) s ( n \ where P n (y 2 ) are polynomials. This results in 
the set of eigenvalues xi n )i 71 > 0, 



X(n) = (tt/2L) 2 <^ 2n + 1/2 + y/W + l/A 



Being substituted into Eq. (f47|), this gives the pre- 
exponential factor J in the main order in Tl r: 

J « exp (£ l= L,fl \ InT, + const) » C (T L T R ) 1/4 . (77) 
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X. APPENDIX A: PRE-EXPONENTIAL FACTOR 
FOR QUASI-1D CASE 

In this Appendix, we show some details of calculations 
of the pre-exponential factor J for quasi-lD case. Due 
to the condition 9 t — at x = x , the observation point 
splits the interval [0, L] into two pieces, and the spectrum 
{x( n )} °f fluctuations around the saddle-point solution 
can be found in each interval separately. Therefore, we 
represent the pre-exponential as a product J — JlJr of 
contributions from the left- and right-hand-side intervals 



XI. APPENDIX B: PRE-EXPONENTIAL 
FACTOR FOR 2D CASE 

In the 2D case, the spectrum {x(n, m)} of fluctuations 
around the saddle-point should be classified by orbital 
and radial quantum numbers n and m, respectively, and 
can be found from the eigenvalues of the Hamiltonian 

H = -r~ x d r {rd r ) + m 2 r~ 2 + U, 

where U is determined by Eq. (fl9|). 

Without any symmetry-breaking, the spectrum of x's 
can be approximated as 

x(0,0)«2L- 2 /ln(L/0 (78) 

for the lowest mode and as 

X (n,m)«(f) 2 [n+i + f] 2 (79) 

for higher n's and m's. 
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The optimal solution breaks the fermion-boson sym- 
metry and induces effective potentials composed of two 
types of contributions, 



-k(d9 t ) 2 fa kfi 2 /r 2 and 



t 



2-nvD 



r 



In the above equations, the approximate values are 
given for the most important range of distances r < 
Ly/z(T)/T, and one has to remember that [i < 1. 

For any m ^ 0, the potential U is relatively small, 
U <C m 2 /r 2 , and could be treated perturbatively. Due 
to the sum rule mentioned in Section III.D, Eq. (p0|), the 
modes with m^O contribute only in the second order in 
U, and what they give to the exponential of J is of the 
order of /j, 4 ln(L/l); /i 2 . With the accuracy we need here 
regarding the leading terms in F t , this contribution can 
be neglected. 

The spectrum of low-lying modes {x( n > 0,0)} is given 
by the expression 



X(n,0) 



0<n<%. (80) 



The cance llation between different eigenvalues from 
substituted to the general equation in Eq. 



Eqs. (78 



( |47| ) produces a multiplier to J which is of the order of 
e~^ 2 for T < 1 and fxhx j at T > 1. Finally, we get 

J = 1 + 0(T 2 ), T< 1 
Jcx /iexp(//ln^) ~ T, T>1" 

This result can be used for all symmetry classes. 



[1] P.W. Anderson, Phys. Rev. B 109, 1492 (1958) 

[2] D. Thouless, Phys. Rev. Lett. 39, 1792 (1972) 

[3] M. Kaveh and N. Mott, J. Phys. C 15, L697 (1982); ibid., 
L707 (1982) 

[4] F. Wegner, Z. Phys. B 35, 207 (1979) 

[5] P.A. Lee and T.V. Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985), and refs. therein 

[6] B. Kramer and A. MacKinnon, Rep. Progr. Phys. 56, 
1469 (1993), and refs. therein 

[7] K.B. Efetov, Adv. Phys. 32, 53 (1983) 

[8] F. Wegner, Z. Phys. B 36, 209 (1980) 

[9] B.L. Altshuler, V.E. Kravtsov and I.V. Lerner, in Meso- 
scopic Phenomena in Solids, ed. B.L. Altshuler et al, El- 
sevier, Amsterdam 1991, p. 449, and refs. therein 
[10] H.Aoki, J. Phys. C 16, L205 (1983); Phys. Rev. B 33, 
7310 (1985) 

[11] M. Schreiber, J. Phys. C 18, 2493, (1985); Phys. Rev. 
B 31, 6146, (1985); M. Schreiber and H. Grussbach, 
Phys. Rev. Lett. 67, 607 (1991); H. Grussbach and M. 
Schreiber, Phys. Rev. B 48, 6650 (1993); ibid. B 51, 663 
(1995) 



[12 
[13 
[14 

[is; 



[16 
[17 

[is; 

[19 

[20 
[21 

[22 

[23 

[24 

[25 

[26 

[27 

[28 



B. Kramer, Y. Ono and T. Ohtsuki, Surf. Sci. 196, 127 
(1988) 

T. Karawarabayashi and T. Ohtsuki, Phys. Rev. B 51, 
10897 (1995) 

S. Evangelou and D. Katsanos, J. Phys. A 26, L1243 
(1993) 

M. Janssen, O. Viehweger, U. Fastenrath and J. Hajdu, 
Introduction to the Theory of the Integer Quantum Hall 
Effect, VCH, Weinheim 1994; W. Pook and M. Janssen, 
Z. Phys. B 82, 295 (1991); U. Fastenrath, M. Janssen 
and W. Pook, Physica A 191, 401 (1992) 

C. Castellani, L. Peliti, J. Phys. A 19, L429 (1986) 
B.A. Muzykantskii and D.E. Khmelnitskii, Phys. Rev. B 
51, 5480 (1995) 

L.P. Gorkov and G.M. Eliashberg, Zh. Eksp. Teor. Fiz. 
48, 1407 (1965) (Sov. Phys. JETP 21, 940 (1965)) 
T.A. Brody, J. Flores, J.B. French, P.A. Mello, A. Pandey 
and S.S.M. Wong, Rev. Mod. Phys. 53, 385 (1981); 
J.J.M. Verbaarschot, H.A. Weidenmuller and M.R. Zirn- 
bauer, Phys. Rep. 129, 367 (1985) 
M. Berry, Proc. Roy. Soc. Lond. A 400, 229 (1985) 
R.A. Jalabert, A.D. Stone and Y. Alhassid, Phys. Rev. 
Lett. 68, 3468 (1992) 

K.B. Efetov and V.N. Prigodin, Phys. Rev. Lett. 70, 1315 
(1993) 

A. D. Mirlin and Y.V. Fyodorov, Phys. Rev. Lett. 72, 526 
(1994) 

B. D. Simons and B. L. Altshuler, Phys. Rev. Lett. 70, 
4063 (1993); Phys. Rev. B 48, 5422 (1993) 

V. N. Prigodin, K. B. Efetov and S. Iida, Phys. Rev. Lett. 
71, 1230 (1993) 

V.I. Fal'ko and K.B. Efetov, Phys. Rev. B 50, 11267 
(1994) 

Y.V. Fyodorov and A.D. Mirlin, JETP Lett. 60, 790 

(1994) 

'Intrinsic' conductance g depends on the dimensionality 
of the system: 

p F l 



ID 



3D : g ~ (p F l)< 



ln(L/0 ' 

Here, we speak about the dimensionless conductance, g = 
G/(e 2 /h). 

[29] W. Feller, Introduction to Probability Theory and Its Ap- 
plications, John Wiley Inc., NY 1966 
[30] B.L. Altshuler and V.N. Prigodin, Zh. Eksp. Teor. Fiz. 

95, 348 (1989) (Sov. Phys. JETP 68, 198 (1989)) 
[31] It is worth mentioning that /(J., t) is normalized, to = 1, 
and the normalization of spin-projected wave functions, 
ti = 1/(2F), produces the same requirement as in Eq. 



[32] V.I. Fal'ko and K.B. Efetov, submitted to Phys. Rev. 
Lett. 

[33] B.B. Mandelbrot, J. Fluid Mech. 62, 331 (1974); T. 
Halsey, M. Jensen, L. Kadanoff, I. Procaccia and B. 
Shraiman, Phys. Rev. A 33, 1141 (1986); G. Paladin and 
A. Vulpiani, Phys. Rep. 156, 147 (1987) 

[34] By this term, we denote the class of the systems both 
with the spin-dependent mechanism of scattering and 
broken time-reversal symmetry. For instance, this can be 
imagined as a metal with built-in disordered static mag- 
netic field or a random exchange field. 



1G 



[35] B. DeWitt, Supermanifolds , Cambridge University Press, 

Cambridge 1992 
[36] V.E. Kravtsov, I.V. Lerner and V.I. Yudson,Phys. Lett. 

A 134, 245 (1989) 
[37] G. Castilla and S. Chakravarty, Phys. Rev. Lett. 71, 384 

(1993) 

[38] F. Wegner, Z. Phys. B 78, 33 (1990) 

[39] K.B. Efetov and A.I. Larkin, Zh. Eksp. Teor. Fiz. 85, 

764 (1983) (Sov. Phys. JETP 58, 444 (1983)) 
[40] A.D. Mirlin and Y.V. Fyodorov, J. Phys. A 26, L551 

(1993) 

[41] R. Feynman, Rev. Mod. Phys. 20, 367 (1948) 
[42] E. Kamke, Differential Equations, Clelsea Publ., NY 
1971 

[43] K.B. Efetov, Zh. Eksp. Teor. Fiz. 92, 638 (1987) (Sov. 
Phys. JETP 65, 360 (1987)) 



17 



